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Abstract 

Using a data sample corresponding to 13.7 fb~ l collected with the CLEO II and II. V detectors, 
we report new branching fraction measurements for two Cabibbo-suppressed decay modes of the 
D+ meson: B(D+ -► vr+vr ) = (1.31 ± 0.17 ± 0.09 ± 0.09) x 10~ 3 and B(D+ -»• K + K s ) = 
(5. 24 ±0.43 ±0.20 ±0.34) x 10~ 3 which are significant improvements over past measurements. The 
errors include statistical and systematical uncertainties as well as the uncertainty in the absolute 
D + branching fraction scale. We also set the first 90% confidence level upper limit on the branching 
fraction of the doubly Cabibbo-suppressed decay mode B{D + — ► K + ir®) < 4.2 x 10 -4 . 



* Submitted to the International Europhysics Conference on High Energy Physics, July 2003, Aachen 
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To lowest order, weak decays of mesons may be described by six quark-diagrams shown in 
Fig. ^ external W-emission, internal W-emission, W-exchange, W-annihilation, vertical W- 
loop, and horizontal W-loop (jlj. When using these diagrams to describe processes, dynamical 
assumptions are often made regarding the relative size of their amplitudes as well as the 
nature of interference terms between diagrams. Measurements of hadronic decays of D + 
mesons give insights into these assumptions as well as new information on flavor SU(3) 
symmetry violation, isospin symmetry, and doubly Cabibbo-suppressed decays. 

Flavor SU(3) symmetry breaking is of current interest because of D° — D° mixing studies; 
it has been shown that the mass and width differences (x, y) of the CP-eigenstates of neutral 
D mesons are generated by the second order SU(3)p symmetry breaking |2j. Understanding 
the size of these effects may be important to unravel any non-Standard Model contributions 
to D°D° mixing. Such understanding is only possible if SU(3)p violating effects are well- 
determined. We report new measurements of the decay modes D + — > 7r + 7r° and D + — > 
K$K + , which are useful for the estimation of SU(3)p violating effects in the D meson 
system. 

Predictions based on isospin symmetry are generally considered to be more reliable than 
SU (3) f predictions because of the near degeneracy in mass of the u and d quarks. Using 
measurements from this analysis as well as data from the Particle Date Group [3], we 
determine the isospin amplitudes and phases for the D — > tttt system. 

Doubly Cabibbo-suppressed decays (DCSD) of charm mesons involve c — > d and s — > u 
quark transitions. Currently, there are only four measured DCSD decay modes 0]. Measure- 
ments of such modes will lead to improved understanding of SU (3) p and other Standard 
Model predictions. Such modes are also important for neutral D-mixing measurements, 
where a significant background is from DCSD decays. In this paper we report the first 
upper limit on the branching fraction of the DCSD decay D + — > K + tt°. 

This analysis uses data collected with two configurations of the CLEO detector at the 
Cornell Electron Storage Ring (CESR): CLEO II Q and CLEO II. V [5]. The total inte- 
grated luminosity of the data sample is 13.7 fb _1 . The CLEO detector is a general purpose 
spectrometer with excellent charged particle and electromagnetic shower energy detection. 
In CLEO II the momenta of charged particles are measured with three concentric drift cham- 
bers between 5 and 90 cm from the e + e~ interaction point. In the CLEO II. V configuration 
the innermost drift chamber was replaced by a 3 layer silicon vertex detector. Charged 
particles are identified by means of specific ionization measurements (dE/dx) in the main 
drift chamber. The tracking system is surrounded by a scintillation time-of-flight system 
and a CsI(Tl) electromagnetic calorimeter. These detectors are located inside a 1.5 T su- 
perconducting solenoid, surrounded by an iron return yoke instrumented with proportional 
tube chambers for muon identification. 

Charged pion and kaon candidates were required to pass minimum track-quality criteria. 
Kaon (pion) candidates had to have a specific ionization within two (three) standard devi- 
ations (a) of that expected for a true kaon (pion). We combined pairs of electromagnetic 
showers in the calorimeter to create 7r° candidates. Candidates with a reconstructed mass 
within 2.5 o of the nominal 7r° mass were kept for further studies. We obtain K$ candidates 
by reconstructing the decay mode K$ — > 7r + 7r~. We required daughter tracks to have an 
impact parameter in the plane transverse to the beam greater than three times the measure- 
ment uncertainty and that the probability of the \ 2 returned from the vertex fit for pairs of 
daughter tracks was required is than 0.001. K$ candidates also had to have a reconstructed 
mass within 3.0 o of the nominal K$ mass. 
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FIG. 1: Six lowest order quark diagrams for a meson decaying into two mesons 0: (a) external 
W-emission, (b) internal W-emission, (c) W-exchange, (d) W-annihilation, (e) horizontal W-loop, 
(f) vertical W-loop. Dashed lines represent W meson. 



In order to reduce backgrounds, we reconstructed the decay process D* + — > D + tt°, re- 
quiring the mass difference (AM) of the reconstructed D* + and D + to be within 2.5 o 
of the known value. We required all D* + candidates to have a normalized momentum 
(xjj* = —j=j=L==) greater than 0.6 and all D + candidates to have a cos9 he [ icity value 

between ±0.8, where OheUdty is the angle between the charged daughter track in the rest 
frame of the D + and the D + in the rest frame of the D* + meson. To insure that we ob- 
tained only one D + candidate per event, we selected candidates with the lowest value for 
y 2 = ( am ~ am pdg) _|_ j2 ( m „o mp. ) w here i indexes the 7r°s candidates in this decay. Given 

A CAM ^ l cr^O ' J 

the large uncertainties in absolute D + branching fractions we present our results as ratios 
of the branching fraction of the decay mode under study to that of a normalization mode: 
D+ -> K-7T+7T+ for D+ -> tt+tt , K + 7i° and D + -> K s tt + for -> 

To extract the yield for each mode, we performed an unbinned maximum likelihood fit 
for two components (signal and background) using the following observables: mo, the mass 
of the reconstructed D + meson, the normalized momentum of the D* + meson, x p m*), and 
cos OheUdty. Using a GEANT-based simulation p of the CLEO detector as well as sideband 
data we determined probability density functions (PDF) for each observable describing the 
shape of the data for signal and background events for each decay mode. The probability that 
a candidate is consistent with signal or background is given by the product of these PDFs. 
The likelihood is given as the product of these probabilities over all candidates; maximization 
of the log of the likelihood gives us the signal and background yields. Projections of the 
likelihood fit to the D + mass for our three decay modes are shown in Fig. |21 Using simulated 
signal and background events, we measure the efficiency of our analysis method for each 
mode, enabling us to determine the total number of signal events in our data sample for 
each decay mode. Table U lists raw yields and efficiencies for all decay modes. 

We considered systematic uncertainties from experimental resolution, efficiency determi- 
nation, and PDF parameterization. The first two contributions are small and the systematic 
errors are dominated by uncertainties in the PDF parameterization. We studied this sys- 
tematic effect for each mode by simultaneously modifying every PDF parameter within its 
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FIG. 2: Invariant mass distributions for D + — ► ir + ir°, K + Ks, K + 7r° cancidates. 



TABLE I: Yields from the maximum likelihood fit with statistical errors and reconstruction effi- 
ciencies. 



Mode 


Yield 


Efficiency 


7T+7T 


171.3 ±22.1 


(6.20 ±0.11)% 


K+K s 


277.7 ± 20.8 


(4.94 ± 0.23)% 


K+TT° 


34.3 ± 20.9 


(6.08 ±0.22)% 




12898.0 ± 156.6 (6.74 ± 0.12)% 




1434.7 ±48.0 


(4.83 ± 0.23)% 



uncertainty. We extracted the yield from the data after each modification to produce a 
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distribution of yields. We defined the systematic uncertainty due to PDF parameterization 
as the 68% limits for these distributions. 

Combining the systematic error study with the yields and efficiencies given in Table U we 
obtain the following results: 

0.0144 ±0.0019 ±0.0010 
0.1892 ±0.0155 ±0.0073 
0.0029 ±0.0018 ±0.0009 



B(D + - 


-> 7T+7T ) 


B(D+ -> 


^-7T + 7T + ) 


B{D+ 




B(D+ 




B(D + - 




B(D+ 


K-71 + 71 + ) 



where the first error is statistical and the second error is systematic. The results supersede 
previous CLEO measurements 0,0. 

In order to determine the absolute branching fractions, we combine our results with the 
PDG values of B{D+ -> K-tx + tx+) = (9.1 ±0.6)% and B(D+ -> n + K°) = (2.77±0.18)% 



and find: 

B(D + -> tt+tt ) = (1.31 ± 0.17 ± 0.09 ± 0.09) x 10" 3 

B(D + ^K + K ) = (5.24 ± 0.43 ± 0.20 ± 0.34) x 10" 3 

£(£>+ -> ir + 7r ) = (2.64 ± 1.64 ±0.82 ±0.17) x 10" 4 

where the third error in the measurements is due to the uncertainty in the normalization 
branching fractions. 

With no significant signal being observed for the doubly Cabibbo-suppressed decay 
D + — > K + 7T° we determined the 90% confidence level upper limit for this branching frac- 
tion. Our method for obtaining the upper limit involved creating 1000 new data sets with 
the same number of signal and background events as our data sample. In order to include 
systematic uncertainties in our upper limit, we also modified the PDF parameters in the 
manner described for our branching fraction calculation. Using this method, our upper limit 
is 

B{D + -> K + 7i°) < 4.2 x 10~ 4 at 90% C.L. 
In the limit of SU(3)p, the following ratio is expected to be unity [7(: 



Rt = 2 x 



Vr, 



cd 



Y(D+ _> tt+ttQ) 
T(D+ -» K°tt+) 



where the V cs and V c d arise because of the different quark transitions in the two decays and 
the factor of 2 arises because of the term in the normaization of the 7r° wavefunction. Us- 
ing the yields and efficiencies (Table Hj) obtained from our analysis and combining statistical 
and systematical uncertainties in quadrature, we find 

Rx = 1.84 ±0.38 

slightly inconsistent with theoretical expectations that SU(3)p symmetry breaking effects 
are about 30%. 
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It is believed that in the D meson system the interference between external and internal 
spectator decay amplitudes is destructive. In order to test this assumption we calculate the 

Ia1l<i 1 T(D+ -f K°K+) _ T(D+ -> K S K+) 

2 ~ 2 X T(D+ 7T + 7T ) ~ r(D+ 7T + 7T ) ' 

which in case of destructive interference should be greater than 1. Besides a small contri- 
bution from the W-annihilation diagram Q the decay in the numerator, D + — > K$K + , can 
be described using an external W-emission diagram. Whereas both the external and the 
internal W-emission amplitudes contribute to the decay in the denominator, D + — > 7r + 7r°. 
Experimentally, we find using our yields and efficiencies from Table U 

R 2 = 2.03 ± 0.32 

indicating that the interference between external and internal W-emission is indeed destruc- 
tive. 

Final state interactions (FSI) are significant in charm decays. Using our measurement 
for D + — > 7T + 7T and the PDG values for D° — > 7r + 7r~, 7r°7r° [3] we can gain some insights on 
these effects by determining isospin amplitudes and phases for the D — > 7T7T system. The tttt 
final state may have an isospin value of or 2. Writing the amplitudes for the 1 = state 
as A and the 1 = 2 state as A 2 , we obtain the following relation: 

A, 2 r +0 



A 



3 



2 (r+- + r 00 ) -r+° 



where T ab = T(D + — ► 7i a TT b ) and a, b represent the charges of the tt. Since isospin amplitudes 
are complex, measuring the phase between them is necessary to obtain full information about 
the amplitudes. The phase is written as 

3r+- - 6r 00 + 2r +0 

coso 



4v / 2T+o^f (r+- + r 00 ) - r+° ' 

We find \A 2 /A Q \ = 0.421 ± 0.044 and cos 5 = 0.042 ± 0.195. These results supersede a 
previous CLEO measurement [8|. The large relative phase between the isospin amplitudes 
indicate that there are significant FSI effects in the D — ► tttt system. A similar observation 
has been made by the FOCUS collaboration (lp| . 

In summary, we have obtained measurements for two singly Cabibbo-suppresed D + decay 
modes: B(D + -> tt+tt ) = (1.31 ± 0.17 ± 0.09 ± 0.09) x 10" 3 and B(D + -> K + K s ) = 
(5.24 ± 0.43 ± 0.20 ± 0.34) x 10~ 3 . We also present an upper limit on the DCSD mode 
B(D + — ► K + 7i°) < 4.2 x 10 -4 at the 90% C.L. Our experimental measurements confirm the 
destructive nature of the interference term between the external and internal W-emission 
diagrams and indicate significant SU(3)p symmetry breaking. An isospin analysis shows 
that FSI effects are important for hadronic decays of D mesons. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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